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Abstract 

Cellular automata are often used to model systems in physics, social sciences, biology that are 
inherently asynchronous. Over the past 20 years, studies have demonstrated that the behavior 
of cellular automata drastically changed under asynchronous updates. Still, the few mathemat- 
ical analyses of asynchronism focus on one-dimensional probabilistic cellular automata, either 
on single examples or on specific classes. As for other classic dynamical systems in physics, 
extending known methods from one- to two-dimensional systems is a long lasting challenging 
problem. 

In this paper, we address the problem of analysing an apparently simple 2D asynchronous 
cellular automaton: 2D Minority where each cell, when fired, updates to the minority state of 
its neighborhood. Our experiments reveal that in spite of its simplicity, the minority rule exhibits 
a quite complex response to asynchronism. By focusing on the fully asynchronous regime, we 
are however able to describe completely the asymptotic behavior of this dynamics as long as the 
initial configuration satisfies some natural constraints. Besides these technical results, we have 
strong reasons to believe that our techniques relying on defining an energy function from the 
transition table of the automaton may be extended to the wider class of threshold automata. 

An abstract version of this paper has been published in /i^ . 



1 Introduction 



In the literature, cellular automata have been both studied as a model of computation presenting 
massive parallelism, and used to model phenomena in physics, social sciences, biology... Cellular 
automata have been mainly studied under synchronous dynamics (at each time step, all the cells 
update simultaneously). But real systems rarely fulfill this assumption and the cell updates rather 
occur in an asynchronous mode often described by stochastic processes. Over the past 20 years, 
many empirical studies [2| [H [5l [131 (EI have been carried out showing that the behavior of a cel- 
lular automaton may widely vary when introducing asynchronism, thus strengthening the need for 
theoretical framework to understand the influence of asynchronism. Still, the few mathematical 
analyses of the eflFects of asynchronism focus on one-dimensional probabilistic cellular automata, 
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either on single examples like [HI [9l [15] or on specific classes like O [7]. As for other classic dy- 
namical systems in physics, such as spin systems or lattice gas, extending known methods from 
one- to two-dimensional systems is a long lasting challenging problem. For example, understand- 
ing how a configuration all-up of spins within a down-oriented external field evolves to the stable 
configuration all-down has only recently been solved mathematically and only for the limit when 
the temperature goes to 0, i.e., when only one transition can occur at time (see [3]). Similarly, 
the resolution of the study of one particular 2D automaton under a given asynchronism regime is 
already a challenge. 

Our contribution. In this paper, we address the problem of understanding the asynchronous 
behavior of an apparently simple 2D stochastic cellular automaton: 2D Minority where each cell, 
when fired, updates to the minority state of its neighborhood. We show experimentally in Section |2] 
that in spite of its simplicity the minority rule exhibits a quite complex response to asynchronism. 
We are however able to show in Section [3] that this dynamics almost surely converges to a stable 
configuration (listed in Proposition [g]) and that if the initial configuration satisfies some natural 
constraints, this convergence occurs in polynomial time (and thus is observable) when only one 
random cell is updated at a time. Our main result (Theorems [s] and [T7| ) rely on extending the 
techniques based on one-dimensional random walks developed in [6l [7] to the study of the two- 
dimensional random walks followed by the boundaries of the main components of the configurations 
under asynchronous updates. We have strong reasons to believe that our techniques relying on 
defining an energy function from the transition table of the automaton may be extended to the 
wider class of threshold automata. 

Our results are of particular interest for modeling regulation network in biology. Indeed, 2D Mi- 
nority cellular automaton represents an extreme simplification of a biological model where the 
biological cells are organized as a 2D grid and where the regulation network involves only two genes 
(the two states) which tend to inhibit each other [Ij. The goal is thus to understand how the con- 
centrations of each gene evolve over time within the biological cells, and in particular, which gene 
ends up dominating the other in each cell, i.e., in which state ends up each cell. Understanding this 
simple rule is thus a key step in the understanding of more complex biological systems. 

2 Experimental results 

This section is voluntarily informal because it presents experimental observations whose formaliza- 
tions are already challenging open questions. The next section will present in a proper theoretical 
framework our progresses in the understanding of these phenomena. The configurations studied 
here consist in a set of cells organized as a n x m torus (n and m are even) in which each cell 
can take two possible states: (white) or 1 (black). The asynchronous behavior of 2D minority 
automaton turns out to be surprisingly complex for both of the studied neighborhoods: 

• von Neumann (N-neighborhood for short), where each selected cell updates to the minority 
state within itself and its neighbors N, S, E, and W; and 

• Moore (M-neighborhood for short), where each selected cell updates to the minority state 
among itself and its 8 closest neighbors N, S, E, W, NE, NW, SE, and SW. 
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Figure 1: 2D Minority under different ^-asynchronous dynamics with N^q = 50 x 50 cehs. The last 
column gives, for a G [0, 1], the empirical probability that an initial random configuration converges 
to a stable configuration before time step Tg • A^5o where Tg = 1000 and Tg = 2000 for von Neumann 
and Moore neighborhood respectively 

In this section, we present a report on extensive experiments conducted on 2D Minority for both 
N- and M-neighborhood. 

In this section, we consider the a- asynchronous 2D Minority dynamics in which at each time step, 
each cell updates to the minority state in its own neighborhood independently with probability 
a. We denote by a = the fully asynchronous 2D Minority dynamics in which at each time 
step, a daemon selects uniformly at random one cell and updates it to the minority state in its 
neighborhood. 

The synchronous regime {a — 1) of 2D Minority has been thoroughly studied in [10] where it 
is proved that it converges to cycles of length 1 or 2. Experimentally, from a random configuration, 
the synchronous dynamics in both neighborhoods converges to sets of large flashing white or black 
regions. 

As soon as a Uttle bit of asynchronism is introduced, the behavior changes drastically for both 
neighborhoods (see Fig. [l] and open our website [12j for animated sequences). Due to the asyn- 
chronism at each step, some random cells do not update and this creates a noise that progressively 
erodes the flashing homogenous large regions that were stable in the synchronous regime. After few 
steps, the configuration seems to converge rapidly to a homogeneous flashing background perturbed 
by random noise. 

Experiments provide evidences that there exists a threshold c^c, ^ -83 and ac ^ .57 for 

the N- and M-neighborhoods respectively, such that if a < o^c, then stable patterns arise {checker- 
hoards and stripes for N- and M-neighborhood respectively). As it may be observed in [12], above 
the threshold, when a > these patterns are unstable, but below and possibly at these 
patterns are sufficiently stable to extend and ultimately cover the whole configuration. 

Convergence in asynchronous regimes. The last column of Fig. [T] shows that experimentally, 
when a ^ the asynchronous dynamics appears to converge at least with constant probability. 
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rapidly to very particular stable configurations tiled by simple patterns known to be stable for the 
dynamics. Above the threshold, when ac < a < 1^ the asynchronous dynamics appears experimen- 
tally to be stuck into randomly evolving configurations in which no structure seems to emerge. 

We will show in Theorem |8] that if at least one of the dimensions is even, the dynamics will almost 
surely reach a stable configuration, for all ^ a < 1, but after at most an exponential number of 
steps. We conjecture that below the threshold ac this convergence occurs in polynomial time on 
expectation if both dimensions are even (the threshold Tg = 2000 is probably too low for the M- 
neighborhood in Fig|T]). We will prove this result in Theorem 17 for the fully asynchronous regime 
under the N-neighborhood under certain natural constraint on the initial configuration. Similar 
results to the ones to be presented below have been obtained in [17J for the M-neighborhood by 
extending of the techniques presented here. 



3 Analysis of fully asynchronous 2D Minority 

We consider now the fully asynchronous dynamics of 2D Minority with von Neumann neighbor- 
hood. Let n and m be two positive integers and T = x the n x m-torus. A n x m- configuration 
c is a function c : T ^ {0, 1} that assigns to each cell {i^j) G T its state Cij G {0, 1} (0 is white and 
1 is black in the figures). We consider here the von Neumann neighborhood: the neighbors of each 
cell (z, j) are the four cells (z ± 1, j) and (z, j zb 1) (indices are computed modulo n and m, we thus 
consider periodic boundary conditions). We denote by = nm, the total number of cells. 

Definition 1 (Stochastic 2D Minority) We consider the following dynamics 6 that associates 
to each configuration c a random configuration d obtained as follows: a cell G T is selected 

uniformly at random and its state is updated to the minority state in its neighborhood (we say that 
cell is firedj; all the other cells remain in their current state: 




1 if Cij + Ci-ij + Ci^ij + Cij-i + Qj+i < 2 
otherwise 



and c^j^i = Cki for all (A:,/) ^ (^, j)- We say that a cell is active if its neighborhood is such that its 
state changes when the cell is fired. 

Definition 2 (Convergence) We denote by the random variable for the configuration obtained 
from a configuration c after t steps of the dynamics: = ^^(c); = c is the initial configuration. 

We say that the dynamics 6 converges almost surely from an initial configuration to a configu- 
ration c if the random variable T = min{t : = c} is finite with probability 1. We say that the 
convergence occurs in polynomial (resp., linear, exponential) time on expectation z/E[r] ^ p(N) for 
some polynomial (resp., linear, exponential) function p . 

As seen in Section |2} any configuration tend to converge under this dynamics towards a stable config- 
uration.^ i.e., towards a configuration where all cells are in the minority state of their neighborhood, 
i.e., inactive. 

Checkerboard patterns. We say that a subset of cells i? C T is connected if R is connected 
for the neighborhood relationship. We say that R is checkerboard-tiled if all adjacent cells in R are 
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Figure 2: Neighborhood's n ame s and transition tables of Minority 6 and its counterpart Outer- 
Totahstic 976 6 (see section 



3.3): only active cells switch their states when fired. 



in opposite states. A horizontal (resp., vertical) band of width it; is a set of cells R = {(^, j) • k < 
i < k -\- w} for some k (resp., R = {(^, j) : k ^ j < k + w}). 



3.1 Energy of a configuration 

The following natural parameters measure the stability of a configuration, i.e., how far the cells of 
the configuration are from the minority state in their neighborhood. Following the seminal work 
of Tarjan in amortized analysis [19j, we define a local potential that measures the amount of local 
unstability in the configuration. We proceed by analogy with the spin systems in statistical physics 
(Ising Model [3J): we assign to each cell a potential equal to the benefit of switching its state; 
this potential is naturally defined as the number of its adjacent cells to which it is opposed {i.e., 
here, the number of cells which are in the same state as itself); summing the potentials over all 
the cells defines the total energy of the configuration at that time. As we consider arbitrary initial 
configuration, the system evolves out- of- equilibrium until it (possibly) reaches a stable configuration, 
thus its energy will vary over time; in particular, as will be seen in Proposition |4j its energy will 
strictly decrease each time an irreversible transition is performed {i.e., each time a cell of potential 



^ 3 is fired). It turns out that this energy function plays a central role in defining, in Section 3.5 
the variant that will be used to prove the convergence of the system. We will see in particular that 
as observed experimentally in Section |2} the system tends to reach configurations of minimal energy 
as one would expect in a real physical system. 

Definition 3 (Energy) The potential Vij of cell {i^j) is the number of its four adjacent cells that 
are in the same state as itself. The energy of a configuration c is defined as the sum of the potentials 
of the cells: E{c) = J2i j ^ij- 

Definition 4 (Borders) We say that there is a border between two neighboring cells if they are in 
the same state, i.e.: 

• the edge between cells {i^j) and {i,j + 1) is an horizontal border if Cij = 

• the edge between cells {i,j) and {i + 1, j) is a vertical border if Cij = 

Definition 5 (Homogeneous regions) An alternating path is a sequence of neighboring cells 
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that does not go through a border, i.e., of alternating states. This defines an equivalence relationship 
« being connected by an alternating path », the equivalence classes of this relationship are called the 
homogenous regions of the configuration. 

Proposition 1 Each homogeneous region is connected and tiled by one of the two checkerboard 
patterns, either ^ or The boundary of each homogeneous region is exactly the set of borders 
touching its cells. 

Proposition 2 The potential of a cell is the number of borders among its sides. The energy of a 
configuration is twice the number of borders. 

A cell is active if and only if at least two of its sides are borders. 

Corollary 3 // both dimensions n and m have the same parity, (Vc) E{c) G 4N; and (Vc) E{c) G 
2 + 4N otherwise. 

The energy of a n x m-configuration belongs to {0, 2,4, .. . ,4A/'} since each pair of adjacent cehs 
in the same state are counted twice and < vij < 4 for ah (i, j). There are two configurations of 
maximum energy 4A^: all-black and all-white. If n and m are even, there are two configurations of 
energy zero: the two checkerboards. If n is even and m is odd, the minimum energy of a configuration 
is 2n and such a configuration consists in a checkerboard pattern wrapped around the odd dimension 
creating a vertical band of width 2 tiled with pattern ffl. 

Energy of stable configurations. A cell is inactive if and only if its potential is ^ 1. It follows 
that the energy of any stable configuration belongs to {0, 2, ... , N}. Stable configurations are thus 
as expected of lower energy If n and m are even and at least one of them is a multiple of 4, there 
are stable configurations of maximum energy N ^ tiled by the "fat"-checkerboard M or ^. 

Energy is non-increasing. Under the fully asynchronous dynamics 5, the energy may not 
increase over time. 

Proposition 4 From any initial configuration c, the random variables E{c^) form a non-increasing 
sequence and E(c^) decreases by at least 4 each time a cell of potential ^ 3 is fired. 

Proof. The variation of the energy of the configuration when the state of a cell of potential v is 
flipped is 8 — 4t' < 0, since active cells have potential ^ 2. □ 

Initial energy drop. Furthermore, after a polynomial number of steps and from any arbitrary 
initial configuration, the energy falls rapidly below 5A^/3, which is observed experimentally through 
the rapid emergence of checkerboard patterns in the very first steps of the evolution: 

Proposition 5 (Initial energy drop) The random variable T = min{t : E{c^) < 5N/3} is al- 
most surely finite and K[T] = 0{N'^). 

Proof. Consider a configuration c with energy E > 5N/3. We will show that either c contains 
a cell of potential > 3 or two adjacent cells of potential 2 in opposite states. Let us proceed 
by contradiction and assume that every cell of c has potential < 2 and that every adjacent cells 
of potential 2 are in the same state. Let 62 (resp. Wi- and W2) be the number of black 
(resp. white) cells of potential ^ 1 and 2. Let consider the bipartite graph that connects each 
black cell of potential 2 to its adjacent white cells of potential ^ 1. Every black cell of potential 
2 is adjacent to exactly 2 white cells of potential ^ 1 and every white cell of potential ^ 1 is 



6 



adjacent to at most 4 black cells of potential 2. The number of edges in the bipartite graph is 
thus at least 262 and at most 4it;i-, it follows that 262 < 4:Wi-. Symmetrically, 2w2 < 46i-. But, 
N = bi- + 62 + '^1- W2 ^ 3(62 + '^2)72, thus the configuration admits at most 2N/3 cells of 
potential 2 and its energy is < + 2A^/3, contradiction. 

Consider now the variant ^^(c^) = 3E{c^)/2 — N^+{c^) where N^+{c^) is the number of cells of 
potential > 3 in For all time t, < ^^(c^) < 6A^. Let N2 be the number of pairs of adjacent 
cells with potential 2 in opposite states. E{c^) is a non-decreasing function of time and each time 
a cell of potential ^ 3 is fired, E{c^) decreases by at least 4; it follows that lK[E{c^^'^) — E{c^)] < 
—4:N^+{c^)/N. A cell of potential 3 may disappear only if itself or one of its four neighbors are fired; 
and each time a cell of potential 2 adjacent to a cell of potential 2 in an opposite state is fired, the 
potential of the later cell increases to 3. It follows that: 

Summing up the two terms yields: 

E[*(c'«) - *(c')l < _^M!a±*W. 

Then, as long as E{c^) ^ 5A^/3, ^^(c^) decreases at each time step by at least on expectation. 
Since is bounded by 6A^, a classic stopping time analysis (see for example. Lemma 2 in [6J) 
shows that after at most 0{N'^) steps on expectation, either E{c^) drops below 5A^/3 or ^^(c^) drops 
below 37V/2 which also imphes that E{c^) < 57V/3. □ 

3.2 Stable configurations 

Proposition 6 (Stable configurations) Stable configurations are the configurations composed of 
checkerboard-tiled bands. More precisely: 

• if n or m is even, the stable configurations are the configurations composed of a juxtaposition 
of horizontal bands (or of vertical bands) of width ^ 2 tiled by checkerboards; 

• if n (resp., m) is odd and m (resp., n) is even, the bands are necessarily horizontal (resp., 
vertical); 

• finally, if n and m are odd, no stable configuration exists. 

Proof. In a stable configuration, every cell touches at most one border. It follows that borders of 
the homogeneous regions form straight lines at least 2 cells apart from each other. □ 

Corollary 7 If n and m are odd, the dynamics 5 never reaches a stable configuration. 

3.3 Coupling with Outer- Totalistic 976 

From now on up to the end of section |3| we assume that n and m are even (with the only exception 
of Corollary [9]). We denote by S the checkerboard configuration of energy defined as follows: 
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a) n and m are even b) only n is even 

Figure 3: Examples of stable configurations. 



= (i + j) mod 2. Given two configurations c and c', we denote by c © the XOR configuration 
c'' such that c^'^- = {cij + c^^) mod 2. 

Dual configurations. As observed above, the fully asynchronous dynamics tends to converge 
from any initial configuration to configurations tiled by large checkerboard regions. It is thus 
convenient to consider instead, the sequence of dual configurations {&) defined by = ^ © c^, in 
which the large checkerboard regions of appear as large homogeneous black or white regions. 
Clearly, the dual sequence & evolves according to the dynamics 5{.) — S©5(S© .), indeed for all t, 
c^+i = ^ © c^+i = ^ © 5{c') = S © 5(S © &) = 5{&). 

By construction, the two dual random sequences (c^) and {&) as well as their corresponding dynamics 
5 and 5 are coupled probabilistically (see [14j): the same random cell is fired in both configurations 
at each time step. A simple calculation shows that the dual dynamics 6 associates to each dual 
configuration c, a dual configuration as follows: select uniformly at random a cell (z, j) (the same 
cell (z, j) as 6 fires on the primal configuration c) and set: 



1 if S > 3 

1 - Cij if H = 2 with H = Ci-ij + q+i j + + 
otherwise 



and c^j^i = Cki for all (A:,/) ^ {i,j)- It turns out that this rule corresponds to the asynchronous 
dynamics of the cellular automaton Outer- Totahstic 976 [TTj. The corresponding transitions are 
given in Fig. [2j 

Stable configurations of Outer- Totalistic 976. We define the energy of the dual configuration 
c and the potentials of each of its cells (i, j) as the corresponding quantities, E{c) and Vij, in the 
primal configuration c. By Proposition [Gj the stable dual configurations under the dual dynamics 
6 are the dual configurations composed of homogeneous black or white bands of widths ^ 2. The 
two dual configurations of minimum energy are all-white and all-black. 

Experimentally, any dual configuration under the fully asynchronous dynamics 6 evolves towards 
large homogeneous black or white regions (corresponding to the checkerboard patterns in the primal 
configuration). Informally, these regions evolve as follows (see Fig. |2|: isolated points tend to 
disappear as well as peninsulas; borders and surrounded points are stable; large regions are eroded 
in a random manner from the corners or bridges that can be flipped reversibly and their boundaries 
follow some kind of 2D random walks until large bands without corners ultimately survive (see 
Fig. lor m)- 
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Figure 4: The coupled evolutions of Minority 6 on the primal configurations (c^) (above) and its 
counterparts Outer- Totahstic 976 (5 on dual configurations (c^) (below). Note that from step 50A^ 
on, (c*) an (c^) are bounded configurations. 

3.4 Convergence from an arbitrary initial configuration 

In this section, we consider arbitrary initial configurations and show that indeed the dynamics 5 
converges to a stable configuration almost surely and after at most an exponential number of steps 
on expectation. 

Theorem 8 From any initial configuration , the dynamics 5 convergences to a stable configuration 
after at most 2N'^^^^ steps on expectation. 

Proof. According to the coupling above, it is equivalent to prove this statement for the dual 
dynamics. The following sequence of (5-updates transforms any dual configuration c into a dual 
stable configuration : 

• Phase I : as long as there are active white cells, choose one of them and switch its state to 
black; 

• Phase II : as long as there are active black cells, choose one of them and switch its state to 
white. 

During phase 1, the black regions expand until they fill their surrounding bands or surrounding 
rectangles. Clearly according to the transition table Fig. [2j after phase 1 of the algorithm, every 
white cell is inactive and thus is either a border or surrounded. In particular, no white band of 
width 1 survived. During phase 2, the black cells enclosed in rectangles or in bands of width 1 are 
eroded progressively and ultimately disapear. Finally, only black bands of width ^ 2 survive at the 
end of phase 2 and the configuration is stable since it is composed of homogeneous white or black 
bands of width ^ 2 (see Proposition [g]). During each phase, at most N cells change their state. 
We conclude that, from any configuration c, there exists a path of length at most 2N to a stable 
configuration. Now, splits the sequence (c^) into segments {c^^^^^ ^ ...^c^^^^^^^) of length 2N . The 
sequence of updates in each of these segments has a probability to be the sequence of at most 

2N updates given above that tranforms configuration c^^^ into a stable configuration. Since these 
events are independent, this occurs after N^^ trials on expectation. We conclude that the dynamics 
5 and thus 5 converge to a stable configuration after at most 2N - N'^^ steps on expectation. □ 

Corollary 9 From any initial n x m- configuration , where n is even and m is odd, the dynamics 
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6 convergences to a stable configuration after at most SN^^^'^ steps on expectation. 

Proof. Consider the cells within the n x (m — 1) rectangle excluding the last column m — 1. Consider 
the dual configuration inside this rectangle and apply the same sequence of updates as above. After 
Phase I, the black regions within the rectangle have been extended up to their bounding rectangles 
and furthermore no proper white horizontal band remains because since m is odd, either one of 
the white cells at the extremity of such a band would be active (whatever the states of the cells in 
the last column are). After Phase II, the black rectangles have been erased as well as the proper 
horizontal black bands (since m is odd, either one of the cells at the extremities of such a band 
would be active). At this stage, the only remaining active cells are within the last column m — 1 
and possibly in either one of the two neighboring columns or m — 2. An extra series of at most 
2n updates allows then to stabihze the cells in these two columns. It follows that a sequence of at 
most 2N + 2n < 3A^ updates stabilize any configuration, which concludes the result by the same 
argument as above. □ 

Example 1 (Conjecture) Draw a rectilinear gray line wrapped twice around the short odd dimen- 
sion of a {2n + 1) X 2n^ -configuration. Cut the configuration along this line and tile the unwraped 
configuration with a checkerboard pattern. Once rewrapped, the only active cells of the configuration 
are along the gray line (see Fig. f^. 



2n+l 




2n^ 

Figure 5: An odd x even configuration with an exponential convergence time? 

The dynamics 5 can converge only after that the gray line is unwrapped, i.e., only after it merges 
with itself somewhere, which is only possible if the line bends itself into a rectangle whose opposite 
corners meet at the same point of the torus. Unfortunately, the ''tension^^ imposed by the wrapping 
around tends to spread apart the two parts of the gray line around the torus (in order to bend itself 
into a rectangle, the n random walks of the corners on the gray line have to synchronize). We 
thus conjecture that this necessary self- crossing of the gray line may only occur after an exponential 
number of steps (which is confirmed by experiments). 



3.5 Convergence from a bounded configuration 

We consider again that n and m are even. We observe experimentally that most of the time, the 
dynamics converges rapidly to one of the two checkerboard configurations of energy zero. We demon- 
strate in this section that if the dynamics reaches a configuration composed of an arbitrary region 
surrounded by a checkerboard, then it will converge to the corresponding checkerboard configuration 
almost surely after a polynomial number of steps on expectation. This corresponds to the analy- 
sis of the last steps of the behavior observed in experimentation. We beheve that the techniques 
developed here may be extended to prove that the dynamics converges to a stable configuration in 
polynomial expected time from any initial configuration (see discussions in section [i]). 
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Definition 6 (Bounded configuration) We say that a configuration c is bounded if there exists 
a (n — 2) X (m — 2) rectangle such that the states in c of the cells outside this rectangle are equal 
to the corresponding states in one of the two checkerboard configurations. W.l.o.g., we assume that 
the upper-left corner of the rectangle is (1, 1) and that the checkerboard is i.e., a configuration c 
is bounded if Cij = (i + j) mod2 for all (i^j) G {(i, j) : (—1 ^ i ^ 0) or (—1 ^ j ^ 0)}. 

Lemma 10 If c is a bounded configuration, S{c) is also bounded. 

Proof. The cells belonging to the checkerboard pattern outside the rectangle have 3 adjacent cells 
in the state opposite to their own states; these cells are thus inactive (whatever the state of their 
adjacent cell inside the rectangle is).n 

A bounded configuration is thus equivalent to a finite perturbation of an infinite planar configuration 
in tiled with the ^ pattern. Since the dual of ^ is the configuration all-white, the dual of a 
bounded configuration is thus equivalent to a finite number of black cells, included into a (n — 2) x 
(m — 2) rectangle within an infinite white planar configuration in 1? . We shall now consider this 
setting. 

Definition 7 (Convexity) We say that a set of cells R Q 1? is convex if for any pair of cells 
and {i + k,j) (resp., + k)) in R, the cells {i + i^j) (resp., {i^j + i)) for ^ £ ^ k belong 
to R. We say that R is an island if R is connected and convex. 

Our proof of the convergence of the dynamics in polynomial time for bounded configurations relies 
on the definition of a variant which decreases on expectation over time. It turns out that in order 
to define the variant, we do not need to consider the exact internal structure of the bounded 
configuration, but only the structure of the convex hull of its black cells. 

Definition 8 (Convex hull of a configuration) For any finite set of cells R G Z^, we denote by 
hull(i?) the convex hull of the cells in R, i.e., hull(i2) = ^\^S C : 5 is convex and S ^ R}. Given 
a bounded dual configuration c, we define the convex hull of c, hull(c); as the dual configuration 
whose black cells are the cells in the convex hull of the black cells of c, i.e., if R = {{hj) • ^ij — 1}; 
hull(c)^j = 1 if and only if (i^j) G hull(i?). We say that a configuration c is convex if c= hull(c). 

We say that c ^ & if for all {i^j), Cij ^ c[j. Let c be a convex dual bounded configuration. We define 
for each black cell (i^j) in c, the island of c that contains cell (i, j); as the maximum connected and 
convex configuration c' such that c'- — 1 and c' ^ c. This defines a unique decomposition into black 
islands of the convex bounded configuration c. 

The variant. We now consider the following variant: $(c) = £'(hull(c))/4 + |hull(c))|, where 
I hu 11(c)) I is the number of black cells in the convex hull configuration hull(c). We will show that 
from any initial configuration c^, $(c*) decreases by at least 1/N on expectation at each time step 
until it reaches the value 0, i.e., until the primal and dual configurations and converge to the 
infinite checkerboard and the infinite all-white configurations respectively. In order to prove that 
$(c^) decreases on expectation, we need to study the evolution of the convex hull of &] for this 
purpose, we introduce a modified coupled dual dynamics 5 that preserves the convexity of a dual 
configuration. Given a dual configuration c, we denote by 5(c) the random configuration c' such 
that: c' = (5(c) if the cell updated by 5 is not a black bridge, and c' = c otherwise. 

Lemma 11 If c is a convex bounded configuration, 5{c) is a convex bounded configuration. 



11 



Proof. The only active transition in 6 that would break the convexity of the black cells is updating 
a black bridge (see Fig. [2|, but this transition is not allowed in 6. □ 

Lemma 12 For all convex bounded configurations c and c' , if c ^ c' , then E{c) ^ E{c'). 

Proof. The energy of a configuration c is by definition twice the number of adjacent cells in opposite 
states in c, that is to say twice the number of sides of cells on the boundaries of the black islands that 
compose c, i.e., twice the sum of their perimeters. Since c < c', the black islands that compose c 
are included within the black islands that compose c^ Moreover, since the sets of rows and columns 
touched by the black islands that compose a convex configuration are pairwise disjoint, the sum of 
the perimeters of the black islands of c that are included in the same black island of c' is bounded 
from above by the perimeter of this later island. □ 

The following lemma proves that the image of the convex hull of c by the dynamics 5 bounds from 
above the convex hull of the image of c by the dynamics 5. 

Lemma 13 For all hounded configuration c, 5{c) < 5(hull(c)). 

Proof. We only need to prove that 1) if 5 updates a white active cell in c, the corresponding cell 
in 5(hull(c)) is black and 2) if 5 updates an active black cell in hull(c), then the corresponding cell 
in (5(c) is white. This is a direct consequence of the coupling of the dynamics of 5 and 5. 

If a white active cell in c is fired and if the corresponding cell in (hull(c)) is white then both cells 
become black. If a white active cell in c is fired and if the corresponding cell in (hull(c)) is black 
then since the cell in c is active it has two black neighbors thus the cell in (hull(c)) has two black 
neighbors. The only kind of active cell with at least two neighbors of the same color under 5 
dynamics is the corner cell. Indeed if a corner white cell in c is black in (hull(c)) then it is a border 
or surrounded cell. Thus if 5 updates a white active cell in c, the corresponding cell in 5(hull(c)) is 
black. 

An active black cell in hull(c) under 5 dynamics is an active black cell in c under 5 dynamics. Thus 
if 5 updates an active black cell in hull(c), then the corresponding cell in (5(c) is white. □ 

Let A$a(c) be the random variable for the variation of the variant after one step of a dynamics A 
from a configuration c, i.e., A$a(c) = $(A(c)) — $(c). 

Corollary 14 For all hounded configuration c, A$^(c) < A$^(hull(c)). 

Proof. By definition, 

A$^(hull(c)) - A$^(c) = (|5(hull(c))| - |hull(5(c))|) + (£;(5(hull(c))) - E {£)))) . 



According to lemma 13, hull(5(c)) < (5(hull(c)) and thus |hull(5(c))| < |(5(hull(c))|. And by Lemma 
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since both configurations are convex, £'(hull((5(c))) ^ £'((5(hull(c))). □ 
Lemma 15 For all hounded configuration c that consists of a unique hlack island^ 

-4/N < E[A$j(c)] ^ -3/A^. 

Proof. Each active cell is fired with probabihty 1/N. According to the dynamics of S (the same 
as the dynamics of S, Fig. [ij except that black bridges are inactive), if c consists of an island of size 



12 



^2, 

E[A$^(c)] = — (#{black corners} + 2#{black peninsulas}) + #{white corners} 
= — -^#{salient angles} + -^#{reflex angles} — —j^-i 



since #{salient angles} — #{reflex angles} = 4 for all convex rectilinear polygon. Finally, if c consists 
of a unique (isolated) black cell, A$j(c) = — □ 

Lemma 16 For any hounded not-all-white configuration c, E[A$^(c)] ^ —ilN, where i is the 
number of islands that compose hull(c). 



Proof. By Corollary 14, E[A$^(c)] ^ E[A$^(hull(c))]. By convexity of hull(c), the sets of rows and 
columns touched by the islands that compose hull(c) are pairwise disjoint. Thus, one can index the 
islands from 1 to ^ from left to right, and the contacts between islands can only occur between two 
consecutive islands at the corners of their surrounding rectangles. Each contact creates at most two 
new active white cells that contribute for each to E[A$j(hull(c))]. The contribution of each 

island to E[A$j(hull(c))] is at most -3/N according to Lemma [iS) It follows that E[A$^(hull(c))] < 
-3£/N + 2{£ -1)/N ^ -£/N. □ 

Theorem 17 The fully asynchronous minority dynamics 6 converges almost surely from any ini- 
tial bounded configuration c to the stable configuration of minimum energy^ and the expected 
convergence time is 0{AN) where A is the area of surrounding rectangle of the black cells in c. 

Proof. Initially and for ah time t > 0, $(c^) < 2(n - 2 + m - 2) + A < 2iV + A. As long as 



& ^ 0, $(c^) > and according to Lemma 16, E[A$^(c^)] < -l/N. It follows that the random 
variable T = min{t : $(c^) < 0} is almost surely finite and E[T] = 0{nA) (by applying for 
example Lemma 2 in [6J); and at that time, and are the stable configurations all-white and ^, 
respectively. □ 

Example 2 (Worst case configurations) Consider the initial dual bounded n x n- configuration 
c consisting of a black 2 x (n — 2) rectangle. The expected time needed to erase one complete line of 
the rectangle is at least Q(nN) = Q(AN). 

Proof. Consider the initial dual bounded n x n-configuration c consisting of a black 2 x (n — 2) 
rectangle. The first time the dynamics 6 will erase a black cell in a given column, this black cell has 
to be a black corner, which was created by the erasure of one of its black neighbors in a adjacent 
column. The expected time between the erasures of the first black cells in a given column and 
of the first black cell in an adjacent columns is thus Q{N) (the expected time to fire the new 
black corner) and the expected time needed to erase one complete fine of the rectangle is at least 
n{nN) = n{AN). □ 



4 Concluding remarks 

This paper proposes an extension to 2D cellular automata of the techniques based on random 
walks developped in [6l[7] to study ID asynchronous elementary cellular automata. Our techniques 
apply as well with some important new ingredients, to the Moore neighborhood where the cell 
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fired updates to the minority state within its height closest neighbors [T7]. We believe that these 
techniques may extend to the wide class of threshold automata, which are of particular interest, 
in neural networks for instance. We are currently investigating refinements of the tools developed 
here, based on the study of the boundaries between arbitrary checkerboard regions in order to try to 
prove that every arbitrary n x m-configuration converges to a stable configuration in a polynomial 
number of steps when n and m are both even (we conjecture a convergence in time 0{N^) for 
non-bounded toric configurations of even dimensions). This result would conclude the study of 
this automaton under fully asynchronous dynamics. The experiments lead in Section [2] exhibit an 
impressive richness of behavior for this yet apparently simple transition rule. An extension of our 
results to arbitrary ^-asynchronous regime is yet a challenging goal, especially if one considers that 
most of the results concerning spin systems or lattice gas (at the equihbrium) apply only to the 
limit when the temperature tends to 0, i.e., when only one transition occurs at a time. 

Acknowledgements. We would like to thank C. Moore, R. D'Souza and J. Crutchfield for their 
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